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ABSTRACT 

We examine the alignment between H2O megamaser disks on sub-pc scales with circumnuclear disks and bars 
on < 500 pc scales observed with HSTfW¥C3. The HST imaging reveals young stars, indicating the presence 
of gas. The megamaser disks are not well aligned with the circumnuclear bars or disks as traced by stars in the 
HST images. We speculate on the implications of the observed misalignments for fueling supermassive black 
holes in gas-rich spiral galaxies. In contrast, we find a strong preference for the rotation axes of the megamaser 
disks to align with radio continuum jets observed on > 50 pc scales, in those galaxies for which radio continuum 
detections are available. Sub-arcsecond observations of molecular gas with ALMA will enable a more complete 
understanding of the interplay between circumnuclear structures. 



1. INTRODUCTION 



Active galactic nuclei have always posed a basic and funda- 
mental problem - how to cram gas that is happily rotating on 



kpc scales onto an accretion disk on AU scales (e.g., Balick & 



|Heckman] ri982). We do not know the mechanism that dissi- 
pates angular momentum and allows gas to accrete. There are 
no shortage of ideas, including major or minor mergers (e.g., 
[Sanders et al.p 988i Ho pkins et al.||2006|), bars or bars w ithin 
bars ( e.g., Shlosman et a l.|1990t|Maciejewski et al."2002: Hunt' 
et al.p OOS; Ki m et al.,2012|), or nuclear spirals (Englmaier & 



Shlosman 2000 i |Maciejewski|2004t [Martini et al.|2003||Ann & 



Thakur|200 5). In a couple of nearby cases, inflows are directly 



observed along circumnuclear spirals in ionized gas on hun- 
dreds of pc scales (e.g., |Storchi-Bergmann et al.|2007[|Davies| 
let al.|2009| ). 

One intriguing clue to the origin of accreting gas comes from 
the well-known misalignment between the rotation axes of ac- 
cretion disks and the galaxy- scale disk, typically measured us- 
ing radio jets as a tracer of the angular momentum on sub-pc 



Space Telescope observations of nine megamaser host galaxies. 
Combining HST/WFC3 imaging from F336W to F160W, we 
are able to identify disk-like structures on < 500 pc scales in the 
majority of the sample galaxies. Even on these circumnuclear 
scales, we find no evidence for alignment with the megamaser 
disks. 

2. CIRCUMNUCLEAR STRUCTURES 

We observed nine megamaser disk galaxies with HST in Cy- 
cle 18 under GO- 12 185. We have examined the M^n-o'^ rela- 
tion of these galaxies already (Greene et al. 2010), and proposed 
to study the correlations between galaxy bulge properties and 
Mbh with HST. The galaxies are all ~ L* spirals withMg ^21 



scales (e.g., |Ulvestad & Wilson, 19 84, ). In this paper we use 22 
GHz water mega maser emission as a tracer of s ub-pc-scale ac- 
cretion disks (e.g.,'Miyo shi et al. 1 1 995 ':'Lo '2005). Early exam- 
ples of megamaser disks revealed the same misalignment, this 
time between the outer galaxy disk and the ^ 0.5 pc molecu- 
lar disk producing the masers (e.g., [Braatz et al.|1997[ |Green 



hill et al. 2009). The interpretation is ambiguous. Perhaps the 



galaxy swallowed some small gas-rich galaxy that is providing 
fuel with a random orientation, or perhaps the accretion process 
involves torques in which angular momentum is not conserved 
as a function of scale. 

The sample of megamaser disk galaxies mapped with Very 
Long Baseline Interferometry has grown considerably in the 
past few years (e.g., Reid et al.|2009t|Kuo et al |2011j , and the 
trend of misalignment continues. Here, we contribute Hubble 
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mag and Hubble types ranging from SO to Sbc ( [Greene et al. 
|2QT51 ). In two orbits, we obtained F336W, F438W, F814W, 
El low, and E160W (roughly UBIJH) images of each galaxy 
with integration times of 1320,430,2140, 150,420 sec respec- 
tively. Two galaxies (NGC 6264 and NGC 6323) are more dis- 
tant than 100 Mpc, which prohibits the analysis presented here. 
Thus we limit our attention to the seven nearest galaxies in our 
sample (Table 1), including five galaxies from Kuo e taL]([2011|), 
NGC 3393 ( [Kondratko et al.|2008| ), and IC 2560 ( .Ishihara et al.| 
|2001| ). Published VLBI maps reveal the orientation of the mas- 
ing disk very precisely in all cases but IC 2560, where we rely 
on unpublished information (Wagner et al. in prep). We fol- 
low the convention of Kuo et al. and refer to the position angle 
(PA) of the megamaser disk as the angle East of North to the 
blue- shifted side of the disk. 

2.1. Identifying circumnuclear structures 

We use a combination of ellipse-fitting and color maps to 
identify organized circumnuclear structures such as nuclear 
bars (e.g., Macieje wski et al . 2002 ; Erwin & Sparke|[2Q03]), 
nuclear rings (e.g., Buta 1986J), or nuclear spirals (e.g.. Mar 



tini et al.||2003| ). The ellipse fits are performed on the E160W 



Table 1. Megamaser Galaxies 
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Note. — Col. (1): Galaxy. Col. (2): Distance to the galaxy (Mpc) from Kuo et al. (2011). Col. (3): Hubble Type, as modified from NED 
based on our HST images. Col. (4): Galaxy diameter (kpc) from NED. Col. (5): Position angle of the galaxy as a whole measured E of N. These 
measurements have an uncertainty of ~ 10°. Our measurements agree with those from NED to within 10° in all cases. Col. (6): Size (kpc) of 
circumnuclear disks or bars. Col. (7): Position angle of structure in Col. (6) E of N. Col. (8): Outer radius of the megamaser disks (pc), as measured 
by Kuo et al. (2012), except for NGC 3393 that was measured by Kondratko et al. (2008). Col. (9): Position angle of megamaser disk measured 
east of north. PA. equals zero when the blueshifted side of the disk plane has zero east offset and positive north offset. Note that all megamaser 
disks except for NGC 4388 have a measured inclination > 86°, and are thus effectively edge-on. All are published but IC 2560, where we rely on 
unpublished information (Wagner et al. in prep). Col. (10): The flux of a detected elongated radio structure (mJy). When there are multiple jet 
structures on different scales, we tabulate the most compact imaged jet. Col. (11): Wavelength of detected jet (cm). Col. (12): The extent of the 
radio jet (pc). Col. (13): Position angle (East of North) of the radio jet (°). Col. (14): Source of the radio observations: (1) Schmitt et al. (2001) 
(2); Ulvestad & Wilson (1984); (3) A. Sun et al. in preparation; (4) Falcke et al. (1998); FIRST: Becker et al. (1994). 



images using the algorithm of |Jedrzejewski| ( |1987| ) as imple- 
mented by the IRAF code ellipse. The code fits ellipses to 
the galaxy isophotes, allowing the center, position angle (PA), 
and ellipticity {e= l-b/a, where a,Z? are the major and minor 
axis respectively) to vary as a function of radius. Circumnuclear 
structures are identified based on peaks and valleys in elliptic- 
ity with corresponding changes in PA pointing to the presence 
of nuclear bars or spirals (e.g., |Erwin & Spa^e||2003| ). Color 
information provides secondary clues as to the nature of the 
structures that we identify. For instance, nuclear rings are iden- 
tified as red or blue rings (e.g., But a & Crocker] 1 993 ). Note that 
while we do not have direct observations of gas yet, blue col- 
ors that identify recent star formation identify gas-rich regions 
by inference. Our primary interest here is not in identifying 
all galactic structures. Instead, we are most interested in the 
projected PA of the circumnuclear (< 500 pc) structures, for 
comparison with the PA of the megamaser disks (Figure [T]). 

We must be careful in interpreting blue colors in the circum- 
nuclear region, because the active galactic nuclei (AGNs) are 
exciting narrow-line region emission that also falls into our 
broad-band filters. We use a pixel-by-pixel fitting method to 
isolate pixels that are dominated in the blue by narrow emis- 
sion lines. For each pixel in the three optical bands (F336W, 
F438W, F814W), we explore a grid of single stellar population 
(SSP) models ( [Bressan et al.||20T2| ), with metallicities in the 
range [Fe/H]= -1 to (although our results are not sensitive to 
the exact range of metallicity), and age running from 1 Myr to 
13 Gyr. Each grid point produces two colors, effectively U-B 
and B-I. Dust is present in almost all of our galaxies, but we 
are only interested in identifying emission-line gas. Thus we 
use, for each pixel and grid point, a range of extinction values 
(Ay = - 50) to modify the SSP colors according to a stan- 
dard extinction law (Girardi et al. 112008 1). We then determine 
the maximum likelihood value for each pixel by comparing the 
likelihoods of all grid points and extinction values. We conser- 
vatively mask all pixels with a low maximum likelihood (corre- 
sponding to x^ > 4). Note that we have not masked the F160W 



images in making our ellipticity profiles. There are emission 
lines in these bands, but their equivalent width is low. As an ex- 
treme example, the famous S -shaped emission region in NGC 
3393 is detected weakly in the F160W image, but has no impact 
on the measured PA, which is nearly orthogonal to the orienta- 
tion of the narrow-line region. 



2.2. Circumnuclear structures and megamaser disks 

In three galaxies (NGC 2273, NGC 2960 or Mrk 1419, and 
UGC 3789) we see evidence for a circumnuclear disk. In the 
case of NGC 2273 and UGC 3789, there is an inner ring which 
manifests as a spike in the ellipticity profile, and then dust-lanes 
interior to the ring most readily identified with spiral arms. In 
the case of NGC 2273, kinematic data confirm this interpreta- 
tion ( Barbosa et al.|2006t|FaTc6n-Barroso et al.|2006| ). We find a 
blue structure in the inner ^3" of NGC 2960 with constant PA 
and ellipticity as the outer disk. Within the inner ^ \" , there is 
a PA change likely associated with the narrow-line region gas, 
as shown (Figure [T]). 

NGC 3393 also shows a bar and maybe a ring on ~ 20^' 
scales, and then has what appears to be a nuclear bar on '^ 3^' 
scales. At the position of the bar, we see an ellipticity maxi- 
mum at a constant position angle, the classic signature of a bar 



(e.g., [Wozniak et~aL]|1995[ [Menendez-Delmestre et al.||2007| ). 
IC 2560 is a difficult case, because the central blue light is 
dominated by narrow-line emission. On kpc-scales (> 30^0. 
the gal axy shows the classic signature of an X-shaped bar (e.g., 
|Li et al.,2011j . However, in the nucleus, dust obscuration and 
narrow-line emission conspire to make identification of any or- 
ganized structure quite difficult. Since the PA profile is flat for 
radii r > 0.5'^ we assign IC 2560 an inner PA of 45°, but cau- 
tion that the identification in this case is uncertain. Finally, we 
have no information for NGC 4388 or NGC 1194, which are 
edge-on, such that dust from the galaxies on large scales ob- 
scures the nucleus. In these two cases we adopt the PA of the 
large-scale disk (Figure |2]). 
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Fig. 1 . — Images of the five megamaser disk galaxies from our HST program that appear to contain circumnuclear disks or bars on < 500 pc scales. North is 
up and East to the left. For each galaxy we show, from left to right, the F814W image on large scales, a three-color image made with F336W, F438W, and F814W, 
and the F160W image. The central and right-most images are shown on the same scale, as indicated by the yellow box in the left-most image. The angular scale 
(in arcsec) is indicated in the upper right of the color image, while the red scale bar in the right-most image indicates 200 pc. In the central three-color image, 
red light indicates both reddened regions and older stellar populations, while blue light indicates recent star formation. The grey regions are pixels that have been 
masked because the colors are dominated by narrow line emission rather than starlight. In the right-most F160W image the red arrows show the orientation of the 
kpc-scale galaxy (G), the megamaser disk (M) and a jet if known (J; see also Table 1). Recall that the megamaser disks are effectively edge-on in all cases, unlike 
the large-scale galaxies, and the arrow points towards the blue-shifted component of the disk. Finally, we show the radial profile from ellipse, along with the radial 
distribution in ellipticity and PA. All identified structures are marked as B=bar, D=disk, N=narrow-line emission, and R=ring. 




Fig. 2. — Three color images for the edge-on galaxies NGC 1194 (left) 
and NGC 4388 (right), constructed as above but without masking. The yellow 
scale bar indicates 1 kpc while the red arrows indicate the jet (J) and the direc- 
tion towards the blue-shifted side of the megamaser disk (M). The blue arrow 
indicates the extended jet (Ext J) directions. 

As an alternate way to look at the data, we also gener- 
ate structure maps from the F814W images ( |Pogge & Martini 
2QQ2| ) shown in Figure [3] Structure maps use deconvolution 



to remove large-scale power and highlight fine structures from 
dust (dark in our maps) and emission-line regions excited either 
by the AGN or star formation (bright in the structure maps). 
These emphasize the dust structures in NGC 2273, NGC 2960, 
and IC 2560. There is dust structure apparent in UGC 3789, al- 
though it is subtle. However, the structure in NGC 3393 is dom- 
inated by the narrow-line region emission (see also the mask in 
Figure [T]). In summary, for the five megamaser disk galaxies 
studied here, whose large-scale galaxies are not edge-on (Fig- 
ure [T]), we find clear evidence for circumnuclear disks or bars 
in four cases. We previously suggested that the majorit y of the 
megamaser disks are found in pseudobulge galaxies (Greene 
et al.|2010| . Now, with the high spatial resolution of HST, we 



have confirmed that there is evidence for ongoing secular evolu- 
tion in the circumnuclear regions of these galaxies, as expected 
for pseudobulges (e.g., ,Kormendy & Kennicuttg2004j . 

Having identified these circumnuclear structures, we now 
must assign a projected position angle to each. We use the PA 
from the ellipse fitting at the scale of the feature, often corre- 
sponding to a local extreme value in ellipticity. As described 
in ,Erwin & Sparke, ( ,2003| ), projection effects (e.g., the super- 
position of a bright bulge and a nuclear bar) can lead to offsets 
between the isophotes and the measured PA. However, based 
on visual inspection of our galaxies, we believe that our values 
are good to ~ 10°. Table 1 contains the orientation and linear 
size of each circumnuclear structure. The nuclear disk in NGC 



2273 (radius of 2" or 260 pc) is well-studied ( Erwin & Sparke 
|2003), and our PA agrees with the literature value within our 
uncertainty of 10°. For the remainder of the paper, we will 
focus on comparisons between the PA of these circumnuclear 
structures on 100-500 pc scales with the megamaser disks on 
sub-pc scales. If we need to refer to larger scales in the galaxy, 
we will refer to the kpc- scale galaxy. 

In principle, we are also interested in the distribution of in- 
clinations for these circumnuclear structures. Taking the ellip- 
ticity measurements at the same radius as the PA measurements 
and assuming thin disks, we find that all of the targets have high 
inclinations (/ > 55°). However, it is not at all clear that we are 
identifying thin disks in all cases. 

We show the distribution of APA=PAdisk-PAmaser in Figure 
|4|2. The first main result of this paper, highlighted in Figure|4|2, 
is that the megamaser disks do not tend to align with the cir- 
cumnuclear structures (^500 pc) traced by HST. We have also 
included NGC 1068, Circinus, and NGC 3079, where the cir- 
cumnuclear disk PA is measured from molecular disks. In the 



case of Circinus, the megamaser disk is warped, and we take 
the PA of the inner disk. The molecular disk in NGC 1068 is 
warped, but in inclination rather than in PA. 

Of course, we are interested in the orientation of the stars 
and the gas. With the HST images, we infer the presence of gas 
based on either blue colors corresponding to recent star forma- 
tion or dust lanes. We have neither gas or stellar kinematic 
observations at these sub-arcsecond scales yet. However, in 
the megamaser galaxies NGC 1068, NGC 3079, Circinus, and 
NGC 2273, detailed gas observations do exist. The detailed in- 
terpretation is different in each case. In NGC 1068, |Schinnerer 
|et al.| pOOO| ) probe ~ 20 pc scales with CO, and show that the 
molecular gas disk has a warp that starts at ~ 70 pc and ex- 
tends all the way to the maser disk on 1 pc scales. NGC 3079 
contains a kpc-scale molecular disk that is aligned with the ma- 
jor axis of the galaxy, and then a 600-pc scale molecular oval 
that is aligned in PA both with the kpc-scale disk and the pc- 
scale masing disk ( |Koda et al . 2002; K ondratko et al.|2005| ). In 
Circinus, the maser disk is nearly perpendicular to the bipolar 
radio jet and a CO outflow ( [Greenhill et al.|2003] ). Finally, the 
kinematic major and minor axes of the disk in NGC 2273 are 
not aligned, and the disk is likely warped (e.g., |Barbosa et aL 
2006| ). The CO disk appears to align with the warped disk on 



200 pc scales f Petitpas & Wilsont2002 ), while the megamaser 
disk is misaligned in PA by 60° . Existing observations of gas 
and kinematics thus suggest that the circumnuclear structures 
observed here are manifestations of torques that facilitate ac- 
cretion onto smaller scales. 

In conclusion, for objects where < 500pc-scale gas or stel- 
lar structure observations are available, 40% show a misalign- 
ment > 20° with the masing disks. We see no tendency for 
alignment between the stellar structures traced by HST and the 
megamaser disks, although the gas structures do tend to align. 
Note that because we focus here on projected position angles, 
we are observing a lower limit on the misalignments. 

3. NUCLEAR JET TO DISK ANGLE 

We have seen that stellar structures on < 500 pc are fre- 
quently misaligned with the sub-pc scale accretion disk. We 
still expect that the nuclear jet should emerge along the rotation 
axis of the sub-pc scale disk (e.g.,I Pringle et al.|1999 ). Since the 
megamaser disks must be close to edge-on in order for ampli- 
fication to occur, the jets should be perpendicular to the mega- 
maser disk positions on the sky. Indeed, Henkel et al. (2005] ) 
found a high rate of megamaser detections in galaxies where 
the jet was likely in the plane of the sky, suggesting that jets are 
usually aligned with the rotation axis of the maser disk. We are 
able to observe the relative orientation of the megamaser disk 
rotation axis and the jet directly. 

Most objects in our sample have been observed in the ra- 
dio continuum at relatively high angular resolution (Table 1), 
includ ing NGC 1194, NGC 3393 ( [Schmitt et al][20QT]), NGC 
2273 pvestad & Wilson' T ^84|), IC 2560 (Morganti et al. 
|T999l ), and NGC 4388 (Falcke et al.|1998j . In the case of UGC 
3789, we simply quote the FIRST fit on kpc scales since no 
higher resolution data are yet available. Finally, for NGC 2960 
we measured a marginally resolved structure at 20 cm from our 
EVLA observations with a size of 20 ± 3^' (A. Sun et al. in 
preparation). Only IC 2560 is unresolved (on the V or 0.2 kpc 
scales measured by Morganti et al.). We caution that these ob- 
servations probe jets on tens of pc scales and in some cases, 
VLBI observations do reveal PA changes on pc scales (e.g.,|Ul- 



N2960 




Fig. 3. — Structure maps (Pogge & Martini 2002) made in with F814W to emphasize the circumnuclear dust structure. Dark regions are dust lanes while bright 
regions correspond to emission line knots due to either star formation or narrow-line region gas. We have zoomed in on comparable regions as seen in the color 
maps in Figure [T] with the black scale bars again indicating 200 pc. In the case of NGC 2273, NGC 2960, and IC 2560, the dust structure is clear. UGC 3789 
appears to show two arms reaching towards the center inside the inner ring. The structure in NGC 3393 is dominated by the narrow-line gas. 



vestad et al.|1998[|Middelberg et al.|2004|). In fact, t he pc-scale 
jet in NGC 3079 is misaligned ("T rotter et"arip^9^ , while the 
diffuse radio emission on pc scales is aligned with the rotation 
axis of the disk ( jDuric & Seaquist|1988| ). Given the alignment 
that we report on tens of pc scales, it would be quite interesting 
to see different behavior on pc scales for more of our sources, 
in light of the possibility that the jet axis may be determined 
on small scales by the black hole spin but on large (tens of 
pc scales) by electromagnetric forces from the disk (McKinney 
|etal.|2013| ). 

In Figure |4|?, we show the distribution of APA for the pro- 
jected jet and megamaser disk angles. Note that the megamaser 
disks are edge-on, so the PA is determined by measuring the an- 
gle on the sky of the observed linear alignment of maser spots 
and is 90° from the rotation axis of the disk. As Figure|4]shows, 
the vast majority of the jets are aligned with the rotation axis of 
the megamaser disk, as expected. We more than double the 



number of systems with a direct disk -jet comparison (.Herm- 
|steinetal.| 19991 IGreenhill et al. 2003 ; [Kondratko et al.l2005| ). 
The largest outHers are are NGC 3079 (APA= 44°) and NGC 
2273 (APA= 60°). In summary, in all but one of the systems, 
the jets on tens of pc to kpc scales are found to align within 
< 15° of the megamaser rotation axis. 

In addition to compact jets, we would eventually like to in- 
vestigate possible kpc-scales jets to seek evidence of a change 
in direction of the accretion disk on longer timescales. For 
instance, radiation-pressure-driven warping may drive sub-pc 
scale disk precession on ~ 10^ year timescales (e.g., Pringle 
1997| ). Furthermore, it is interesting to ask whether ioniza- 



tion cones align with the jet/megamaser disk rotation axis. We 
do not really have complete information on the orientation of 
the narrow-line regions for the full sample. However, from the 
masked (grey) pixels in Figure [T] we can see that the spectac- 
ular S-shaped narrow-line region in NGC 3393 does align with 
the radio jet and the megamaser rotation axis (e.g., |Cooke et al.| 
|2000t [Kondratko et al.|2008j , and this is seen in Circinus as 
well ( [Greenhill et al.|20Q3 ). In other cases, we do not have ad- 
equate information yet to address the relative orientation of any 
ionization cones. 

4. WARPS AND KINEMATIC MISALIGNMENTS FACILITATE 
ACCRETION 



We have seen that nuclear (and in some cases extended) ra- 
dio jets align with the rotation axis of the megamaser disks. 
In contrast, the sub-pc scale accretion disk does not align with 
the kpc-scale galactic disk (see Q. Using HST /WFC3 imag- 
ing, we have found no strong tendency for alignment between 
the megamaser disks and the circumnuclear structures traced by 
stars on < 500 pc scales. Here we review what these misalign- 
ments may tell us about the accretion process. 

Disk warping on pc scales causes the observed misalign- 
ments. One route to disk misalignment is warping of the pc- 
scale accretion disk. According to Pringle (e.g., 1997), the sub- 
pc scale disk may change position angle by a significant fraction 
due to radiation pressure. For the megamaser disks, the typical 
timescale for disk PA changes due to radiation-pressure warp- 



ing would be ^ 5 x 10^ yr (see also Maloney et al.|l996| 



Gam- 



mie et al.]|2000| ) with only a linear dep endence on BH mass 



AGN lifetimes are very uncertain, but [Martini et al.| ( |2003| ) 
suggest that accretion episodes in low-luminosity sources may 
only last millions of years, perhaps just long enough for some 
PA change over the AGN lifetime. Note that radiation pres- 
sure warping cannot cause the warps observed on circumnu- 
clear (- 200 pc) scales in NGC 2273 and NGC 1068. 

Another possible source of sub-pc-scale torques comes from 
stars in a cusp around the BH; [Bregman & Alexander] ( |2012| ) 
find that resonant relaxation with a stellar cusp will drive warp- 
ing of order 10° on sub-pc scales, comparable to the warp in 
NGC 4258. The expected timescale for evolution is ~ 10^ yr 
for NGC 4258, perhaps in some tension with the million-year 
lifetimes from Martini et al. Alternatively, JKartje et al.| ( [l999| ) 
suggest that masing clumps are accelerated out of the disk plane 
by magnetic pressure to a height where the balance of shielding 
and pumping is optimal. Since the optimal height will depend 
on the AGN luminosity, they predict more warping at higher lu- 
minosities. Both radiation pressure and stellar resonance warp- 
ing act on sub-pc scales, i.e. on the scales of the megamaser 
disks. Warping does not appear to be a compelling explana- 
tion in general for the observed misalignments for three rea- 
sons. First, in contrast to NGC 4258, t he observ ed warps in the 
new megamaser disks are < 10° (Impellizzeri et al. 2012} Reid 
|et al.||2013t |Kuo et al.|r2013i ) and thus smaller than predicted 
by either warping model. Second, the observed warps are too 
small to explain the broad range of PA differences that we ob- 
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Fig. 4. — Histogram of the difference in projected position angle between the megamaser disk and (left): the circumnuclear structures on <500 pc scales identified 
in the HST images and (right): the most compact jet structure we could find in the literature, as summarized in Table 1 (note - this is the PA of the megamaser 
disk, 90° from the rotation axis). Left: In the case of NGC 1 194 and NGC 4388 we simply use the large-scale PA of the galaxy. The filled histograms indicate the 
literature megamasers, where we use the difference in projected position angle between known CO disks in NGC 1068 ( Schinnerer et al. 2000), Circinus (Curranj 
|etal.|199 8>, and NGC 3079 ( Kondratko et al. 2005 1. Thus, these measurements are not strictly comparable (ours probes recent star formation while the other directly 
probes molecular gas). Right: In this case the filled histograms include NGC 4258 as well. 



serve. Third, the timescales for these warps to operate may be 
long compared with typical lifetimes. 

Changes in angular momentum as a function of scale result 
naturally in accretion events. Hopkins & Quataert| ( |2010| ) study 
the progression of gas from galaxy-wide to pc scales. Because 
the gas is dissipational, the kinematics of the gas and stars de- 
couple. As a result, the stars are able to torque the gas, leading 
to warps and ultimately shifts in the angular momentum vector 
of the gas as a function of scale. Observations of the gas kine- 
matics on ~ 100 pc scales in nearby active galaxies reveal that 
fueling is indeed facilitated by decoupled dynamical compo- 
nents such as inner bars, ovals or spirals (e.g.. Hunt et al.|2008[ 
Garcia-Burillo et al. 2009), but that these fueling ep isodes are 
stochastic and short-lived (e.g., Dumas et a l. 2007; Haan et al. 
2009| ), in keeping with the observed kinematic misalignment on 



all spatial scales that we can probe. 

In the Hopkins & Quataert models a massive bulge compo- 
nent will suppress these torques, thus requiring different feed- 
ing mechanisms (e.g., merging) in bulge-dominated system s 
(see also their analytic model in Hopkins & Quataert 2011 ). 
It is intriguing to note that the megamaser aligns with the large- 
scale disk in the two bulge-dominated galaxies NGC 1194 and 
NGC 2960. Our statistics are obviously quite limited. We can 
use our strong confirmation that nuclear jets align with the ro- 
tation axis of megamaser disks (and by extension sub-pc scale 
disks in general). We look for morphology-dependent align- 
ment in the |Kinney et al.|pOOO| ) Seyfert sample, with morpho- 
logical types from the RC3. There is no preference for jet align- 
ment with the rotation axis of the kpc- scale galaxy disk among 
SO or so/a galaxies. However, it would be interesting to inves- 
tigate possible jet alignment with the rotation axes of circum- 
nuclear ~ 500 pc-scale structures in the Kinney et al. sample 
(e.g., using HST). 

Disk misalignments boost accretion rate. I Nixon et"ar ( 2012 ) 
show that if there are strong misalignments between inner and 
outer disks, then the gas angular momentum is dissipated where 
the disks meet, facilitating accretion. As a result, if there is a 
major misalignment as a function of scale, then the accretion 
rate onto the BH can be boosted by an order of magnitude. 



We may preferentially observe AGN activity at times of mis- 
alignment. With a much larger sample, we could look for an 
anti-correlation between AGN luminosity and PA alignment as 
predicted by this model. 

The gas has an external origin, and thus knows nothing about 
the angular momentum of the disk. It would be very challeng- 
ing to completely rule out that the gas is supplied at random 
angles via accretion of small satellites. In fact, in the case of 
NGC 3393 there is tentative evidence for an accretion event in 
the putative detection of two AGN separated by ~ 150 pc in 
projection ( IFabbiano et al.|20lT] ). Nevertheless, none of the tar- 
gets shows evidence for morphological disturbance. Given the 
high incidence of circumnuclear disks and bars in the sample, 
an internal origin for the accreted gas, with accretion facilitated 
by secular processes, seems more natural for the bulk of the 
sample (e.g., |Kormendy & Kennicutt|2004 ). 



5. SUMMARY 

Using HSTfWFC3 observations of megamaser disk galaxies, 
we have shown that the majority of the megamaser disk galax- 
ies contain clear evidence for circumnuclear structures (disks 
or bars) often with associated star formation. We find that the 
accretion disk on sub-pc scales (as traced by the megamaser 
disks) shows no strong tendency to align with circumnuclear 
(often star- forming) structures such as bars, spirals, or rings on 
< 500 pc scales. In contrast, we find that the rotation axis of 
the megamaser disk is usually aligned with the compact radio 
jet axis. We review possible explanations for the misalignment 
between circumnuclear and sub-pc scale disks, and favor a sce- 
nario in which accretion naturally requires changes in angu- 
lar momentum as a function of scale. Given the good correla- 
tion between the jet and megamaser rotation axis, in the future 
we can boost the sample statistics using HST observations of 
Seyfert galaxies with radio jets. Furthermore, future direct ob- 
servations of circumnuclear gas in the megamaser disk galaxies 
(e.g., with ALMA) will help clarify the processes that feed the 
central monster. 

We would like to thank Paul Martini, Jeremy Goodman, Eliot 



Quataert, Jim Ulvestad, and John Kormendy for very useful 
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